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The selection of the ground state among nearly degenerate states due to quantum fluc- 
tuations is studied for the S = 1/2 XY-like Heisenberg antiferromagnets on the triangular 
lattice in the magnetic field applied along the hard axis, which was first pointed out by 
Nikuni and Shiba. We find that the selected ground state sensitively depends on the degree 
of the anisotropy and the magnitude of the magnetic field. This dependence is similar to 
that in the corresponding classical model at finite temperatures where various types of field 
induced phases appear due to the entropy effect. It is also found that the similarity of the 
selected states in the classical and quantum models are not the case in a two-leg ladder 
lattice, although the lattice consists of triangles locally and the ground state of this lattice 
in the classical case is the same as that of the triangular lattice. 

KEYWORDS: Frustration, Triangular lattice, Entropy-induced phase, Quantum fluctuation in- 
duced phase, PWFRG, DMRG 

1. Introduction 

The ordering process in triangular lattice antiferromagnets has attracted much interest 
because the frustration among antiferromagnetic interactions causes complicated order pa- 
rameters. In such systems, various spin configurations are nearly degenerate. Among them 
one configuration is selected as a result of a very subtle balance of the energy and the entropy. 
At finite temperatures, a structure that is entropically favorable would be chosen even if it 
is less favorable in energy. In such cases, we have a phase transition between the entropically 
favorable ordered phase and an energetically favorable ordered phase. Indeed, such fluctuation- 
induced phase transitions have been studied in the magnetization process of antiferromagnets 
on the triangular lattice.^"'' 

A similar selection of the states occurs also for the ground state of a quantum triangu- 
lar lattice antiferromagnets at zero temperature; the classical ground state is not necessarily 
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most stable and some other configuration may become tire ground state wlicn a contribution of 
quantum fluctuations to the energy is taken into account. Among nearly degenerate classical 
configurations, a certain state is selected, where quantum fluctuations play a role similar to 
thermal fluctuations in the classical case at finite temperatures. In particular, such example 
was pointed out by Nikuni and Shiba^"^-*^ to understand the ground state phase transition of 
CsCuCls.^^' ""^^ This material has a hexagonal lattice consisting of layers of triangular antifer- 
romagnets with a weak XY anisotropy and ferromagnetic couplings between the layers. It has 
also been pointed out that RbFeCla has qualitatively the same structure, and a similar phase 
diagram has been obtained. 

In the classical spin model, the ground state of the XY-like Heisenberg antiferromagnet on 
the triangular lattice consists of the 120° structure. This structure of the ground state causes 
the so-called chiral phase transition due to its two- fold discrete degeneracy. When the field 
is applied in the z-direction, the 120° structure is folded to be a shape of umbrella which is the 
most energetically favorable in the field. (Fig. 1(a)) The magnetization increases proportionally 
to the field strength up to the saturation. There is no singularity in the magnetization process. 
However, a jump of the magnetization was found in the experiment of CsCuCls. Nikuni and 
Shiba have pointed out that this jump is due to a phase transition in the ground state. The 
quantum correction on the ground state energy can cause a ground state phase transition. 
By comparing energies of the umbrella structure and a coplanar structure depicted in (Fig. 
1(b)) which we shall call 'v-shape' structure,*^' ""^^ they found that the umbrella structure is 
favorable at low field while the v-shape structure is favorable at high field. That is, although 
the latter has higher energy in the classical system, it has the lowest energy if the quantum 
correction of the energy is taken into account. 

Recently, Watarai et al.^ have pointed out that the role of quantum fluctuations in the 
above scenario for the quantum case is very similar to that of thermal fluctuations in the 
corresponding classical spin system at finite temperatures. Moreover they have shown on the 
basis of the detailed temperature-field phase diagram of the classical system that the structure 
of the phase diagrams is sensitive to the degree of the anisotropy. That is, no phase transition 
occurs when the XY anisotropy is strong, and the phase transition from the umbrella structure 
to the v-shape phase occurs when the system has a moderate anisotropy. In the Heisenberg 
limit, the system shows a qualitatively different sequence of phases from that of the case 
of the moderate XY anisotropy. In the Heisenberg system, the ground state is continuously 
degenerate.^ In the low-field region a coplanar structure, which we shall call a 'Y-shape' 
structure (Fig. 1(c)), is selected by the entropy effect. In the high-field region, the phase of 
the v-shape structure becomes the most stable state as in the case of the XY-like Heisenberg 
model. Between the low and high field phases, a collinear phase (Fig. 1(d)) is present, which 
leads to the 1/3 plateau in the magnetization process. When the XY anisotropy becomes very 
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weak, a complex sequence of phases appears including phases in both of the above two cases. 
Then, it is an interesting question how such sensitive anisotropy dependence in the classical 
model can reflect on the corresponding quantum spin system with various XY anisotropy. 

Although from a view point of the spin wave theory, quantum fluctuations should cause 
an effect similar to that of thermal fluctuations,^*^ it is still interesting to sec how it appears in 
strongly fluctuating quantum systems. In this paper, we study whether similar phase structure 
with strong dependence on the anisotropy appears or not by investigation of the ground 
state properties in the magnetic field by various numerical methods. In the case of three-leg 
ladder system, which well represents the triangular lattice, the dependence is certainly similar 
to those found in the corresponding classical systems at finite temperatures as the entropy 
effect. We also study the dependence of the ground state on the lattice shape. In the classical 
system, the ground state configuration of all the lattices consisting of local triangles is an 
assembly of local ground state configuration of three spins representing the three sublattices 
of the triangular lattice. Thus, the ground state configurations of those lattices have the same 
spin configuration. By contrast, the interaction in quantum systems is not local because of 
the non-conimutativity among the spins, and the ground state configuration depends on the 
lattice shape even if they locally consist of the same triangles. In fact, we demonstrate such 
dependence in a two-leg ladder system. 

In order to study the quantum ground state of the present model, we cannot apply the 
quantum Monte Carlo simulation because of the negative-sign problem. Therefore, in this 
paper, we adopt kinds of diagonalization method, i.e., exact diagonalization method for a finite 
lattice of small number of spins, the product wave-function rcnornialization group (PWFRG) 
method which is a variant of density matrix rcnormalization group (DMRG) method. 
PWFRG has been developed to study the magnetization process in the ground state. We 
also use the DMRG method for finite lattices which gives complementary data for those from 
PWFRG. 

This paper is organized as follows. We explain the model studied in the present paper and 
order parameters in §2, where corresponding data in the classical model are also reviewed. In 
§3, the dependence of the ground state configuration on the anisotropy and on the external field 
is studied in a finite lattice. In §4, the dependence is studied in a two-leg ladder system which 
is an assembly of triangles. In §5, properties on a three-leg ladder with the same periodicity 
as the triangular lattice is studied. In §6, summary and discussion are given. 

2. Model and Order parameters 

For the study of CsCuCl3, three-dimensional XY-like Heisenberg model has been adopted. 
In the classical model, the anisotropy is usually introduced by the single-ion anisotropy 

H = Y. -J^^^^ ■ - ^Y^^Stf -H^S!, (2.1) 

(ij) i i 
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(a) 



(b) 



(c) 



(d) 



Fig. 1. Candidates for the ground state configurations: (a) umbrella structure, (b) v-shape structure, 
(c) Y-shape structure, and (d) coUinear ferrimagnetic structure. 



where {ij) denotes the interaction in a triangular layer. The interaction in the layer is antiferro- 
magnetic, i.e., Jij = J > 0, while the interaction for the interlayer interaction is ferromagnetic, 
i.e., Jij = J' < 0. For simplicity we take J' = — J. 

In the case of = 1/2 quantum system, however, the single-ion anisotropy does not exist. 
Therefore, we adopt anisotropic coupling as a source of the anisotropy. In this paper, we study 
only two-dimensional systems because of the limit of the capacity of numerical calculation. 
The Hamiltonian is given by 

n = J2 J^j {sfSj + s^sy + Asfsp^ -hY^si, (2.2) 

{ij) « 

where Sf = \(yf, {a = x,y and z) where af is the a component of the Pauli matrix. The 
magnetization of the system is measured by 

M = ^5f. (2.3) 

i 

Low-energy configurations of the model (1) in the classical case is depicted in Figs.l, where 
the angles 6 are respectively, given by 

cos^ =3(2J+T) for Fig. 1(a), (2.4) 

cos^ for Fig. 1(c). (2.5) 

For the configuration Fig. 1(b), we have to calculate 6 numerically. Comparing energies of the 
configurations, we find that the umbrella configuration (Fig. 1(a)) gives the minimum energy 
for A < 1, corresponding to D < in the classical spin model (2.1). 

In order to characterize the phases, we define the following order parameters. The unit 
cell in the ordered phases consists of three sublattices. First, we introduce the chirality:^^ 

K=-^ 5Z (SixSj + SjxSk + SkxSi), (2.6) 

where i,j and k denote sites on a triangle consisting of the sublattices A, B, and C of the 
triangular lattice, respectively. The presence of the umbrella structure can be identified by 
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measuring the z component of the chirahty. 

In the v-shape structure, the spins on the two of sublattices are parahel, and thus the chi- 
rahty disappears. In order to identify noncolhnear ordering structures of the v-shape structure, 
wc measure the fohowing quantities: 

ieAjeB ieBjeC igCJgA 

We define the following order parameter to detect the non-collinearity 

Xc^w^ [i^^y + (x^y + ix^^] , (2.8) 

where a = x,y, and z, and N is the number of spins. In particular, we use an order parameter 



Sxy — Xx + Xy (2-9) 

to detect the v-shape structure when the field is applied in the z direction. 
If the configuration is of the umbrella type, the order parameters 

K, = ^{kI) (2.10) 

and Sxy are nonzero in the limit N ^ oo. On the other hand, in the v-shape structure, is 
zero while Sxy remains non-zero. 

The xy component of the chirality K is good to detect the Y-shape structure. Therefore, 
we also observe the following quantity: 

K.y^^{{4) + {4))- (2-11) 
This order parameter is non-zero in the Y-shape structure, but it is zero in the v-shape 
structure. 

For the classical model (2.1), we can perform Monte Carlo simulations in large systems 
and determine the thermodynamic stable state. In Figs. 2, the field-dependence of the chirality 
and Sxy of a weakly XY-like Heisenberg model, i.e., D/J = —0.1, at temperature T/J=0.5, 
are depicted. There, as the field increases, the structure changes from the umbrella to the 
v-shape discontinuously and finally to the field-induced ferromagnetic state. The data for 
the Heisenberg model are also depicted in Figs. 2(c) and (d).^ There the structure changes 
from the Y-shape, collinear up-up-down, and the v-shape, and finally to the field-induced 
ferromagnetic state. 



3. Observation on a finite lattice 

For the present quantum systems, Monte Carlo simulation is not available because of 
the negative-sign problem. Thus, the most naive approach to study the ground state is a 
diagonalization of the Hamiltonian of a finite lattice of small number of spins. In this section. 
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(a) 
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(c) (d) 

Fig. 2. The field dependence of (a) Kz and (b) S^y at T = 0.5 for the weakly XY-anisotropic case 
{D = —0.1). (c) Kxy and (d) Sxy at T = 0.5 for the Hciscnbcrg model (D = 0). The symbols, 
circle, square, diamond, and (closed) triangle denote the size L = 12, 18, 24, and 30, respectively, 
(from the reference^) 



we study systems in a lattice of 3 x 6 with the periodic boundary condition depicted in Fig. 3(a) . 
For the finite size lattice, we study the property of the state with a fixed magnetization M 
instead of the magnetic field H. The change of M corresponds to that of H because the 
magnetization of the ground state monotonically increases with the field. We note that, in 
Refs.,^^'^^ magnetization processes of the similar systems have been also studied for finite size 
and/or ladder lattices. 

In a finite lattice, we investigate correlation functions of the quantities corresponding 
to the order parameters, because the expectation values of the order parameters themselves 
may vanish in the quantum ground state. Here we observe the correlation functions of the 
z-component of the chirality (CKz), the xy-component of the chirality (CKxy), and the xy- 
component of the vector product of spins defined in (2.7) (CSxy), at the sites denoted in 
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Fig. 3. (a)The finite lattice for the diagonalization study. The correlation functions of the order 
parameters are obtained at the quantities on the positions denoted by crosses. The field dependence 
of (b) CKz and (c) CK^y, and (d) CS^y as functions of the magnetization M. Symbols denote 
the data for various values of A: •(A=0.0),o(0.1),®(0.2),(g)(0.3), 0(0.4), <(0.5),v(0.6), A(0.7), 
0(0.8), □(0.9), and O(l-O) 

Fig. 3(a). CKz,CKxy, and CSxy are given by 

CK^ = {[Si X S2 + S2 X S3 + S3 X Si]^ • [S^/ X +S2' X S3/ +S3/ X S^/] J (3.12) 
CKxy = {[Si X S2 + S2 X S3 + S3 X Si]^ • [S^/ X +S2' X S3/ +S3/ X S^/]^) 

+ ^[Si X S2 + S2 X S3 + S3 X Si]y- [S^/ X S2' +S2' X S3/ +S3/ X S^/]^^ (3.13) 

CSxy^{[Si X S2],-[S,, X 82,]^.+ [S2 X S3],.-[S2. X 83,]^,+ [S3 X Si]^,[S3, X S,,] J 

+ ([Si X S2],-[S,, X 82']^+ [S2 X S3],-[S2' X 83,]^+ [S3 X Si]^-[S3, X S,,]^),(3.14) 

where 1,2,3 and 1 ,2 ,3 are the lattice sites of the triangles denoted by crosses. We plot 
CKz, CKxy and CSxy as functions of the magnetization for various values of the anisotropic 
coupling A [A = 0,0.1, • • • , 1) in Fig. 3(b), (c), and (d), respectively. 



7/15 



J. Phys. Soc. Jpn. 



Full Paper 




V 



+ 



A 



Fig. 4. Two-leg ladder consisting of triangles. The correlation functions of the order parameters are 
obtained at the quantities on the positions denoted by crosses. 

We find a qualitative similar dependence of the order parameters to those in the classical 
models. That is, CK^ decreases with the magnetization M smoothly for the case of strongly 
anisotropy {A = 0), where no transition occurs until the saturation field Hco- 

When the system has weak anisotropy, e.g. A = 0.9, then CK^ decreases suddenly at an 
intermediate value of the magnetization M, although CSxy is still large. This case corresponds 
to the classical case depicted in Figs. 2(a) and Figs. 2(b). At the Heisenberg case {A = 1), 
instead of CK^, CK^y appears at small magnetizations, and it disappears at an intermediate 
field where the collinear up-up-down state appears. For larger M, another non-collinear state 
appears. These dependences agree with those in Figs. 2(c) and Figs. 2(d). Thus, the anisotropy 
and field dependence of the ground state property is similar to that in the classical systems, 
although number of data points are limited because only the finite lattice is studied. In the 
next section, we study systems consisting of infinite number of triangles. 

4. Observation on a ladder system consisting of triangles 

In this section, we study a two-leg ladder system with zigzag structure, which consists of 
triangles depicted in Fig. 4. The Hamiltonian of this case is given by 



Heisenberg antiferromagnets with A = 1 in this type of lattice have been studied in detail 
as the Majumdar-Ghosh model. There, various interesting properties, e.g., energy gap be- 
tween the ground state and the first excited state, and cusp singularities in the magnetic 
process, have been found. Here, we study properties of the XY-like antiferromagnets on this 
lattice {A < 1). For this type of one dimensional lattices, methods of the DMRG type are 
very efficient. Here, we study the magnetization process by PWFRG and DMRG. As in the 
exact diagonalization for a finite size lattice, the expectation values of the order parameters 
themselves may vanish in the quantum case. We, therefore, measure the correlation functions 
CKz and CSxy between two triangles near the center of the ladder, which are indicated by x 
symbols in Fig.4. 





(4.15) 
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Fig. 5. The field dependence of the magnetization m, CK^ and CSxy as functions of the field H for 
(a) A = 0.9, (b) A = 1.0 and (c) A = 0.0 



In Fig. 5(a), we depict the magnetization for A = 0.9 obtained by PWFRG. In the present 
study we use a block Hamiltonian of the dimension m = 48. We have confirmed that the results 
converge up to this value of ra. In the magnetization process, we find a zero magnetization 
region at low field, which implies the existence of energy gap. This gap is due to the dimering 
effect that brings the Majumdar-Ghosh solution in the Heisenberg limit. In Fig. 5(a), the 
correlation functions of the order parameters CKz and CSxy are also plotted. They behave 
very difi^crcntly from those expected in the classical lattice. This may be due to the fact that 
the zigzag chain does not have the three sublatticc structure in the rung direction. For example 
CKz takes a negative value for ^ > 0, for which there is no classical analogue. Thus, in the 
next section we study a lattice which is more close to the triangular lattice. 

Before proceeding to the next section, here we make a comparison between the results for 
PWFRG and DMRG. In PWFRG, we obtain the magnitude of the magnetization for the given 
field. On the other hand, in the DMRG procedure, we obtain a ground state energy for given 
values of the magnetization from which we can obtain the magnetization curve as a function 
of the field because the magnetization commutes with the Hamiltonian for the present model. 
Both methods should give the same results. In Figs. 6, we demonstrate that the PWFRG 
and DMRG results agree with each other by showing the quantities for the Majumdar-Ghosh 
model of L = 128 system obtained by the DMRG with m = 64. In the figure we taJie spatial 
averages of CK^ and CSxy for the comparison with the PWFRG results, since CKz and CSxy 
in a finite size system show position depending behavior due to the open boundary. 

5. Three-leg triangular lattice 

In this section, we perform PWFRG on a 3 x oo lattice with the periodic boundary con- 
dition depicted in Fig. 7. We expect that the model on this lattice reproduces well properties 
of the model on the triangular lattice. Since this lattice contains the minimum length of the 
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(a) (b) (c) 

Fig. 6. The field dependence of (a) the magnetization m, (b) CK^ and (c) CSxy as functions of the 
field at A = 1 obtained by PWFRG and DMRG. 




Fig. 7. The 3-leg ladder. The correlation functions of the order parameters are obtained at the quan- 
tities on the positions denoted by crosses. 



unit cell of the triangular lattice both in the rung and leg directions. The unit cell of the order 
parameter is 3 x 3, and therefore it had better take this unit cell as a unit of the renormaliza- 
tion process. However, it is difficult to perform this process because of the limitation of the 
memory. Thus, we take a 3 x 1 cell as a unit cell of the RG process. In Fig. 8, the procedure 
of PWFRG is illustrated. When we renormalize the block Hamiltonian by one step, the sub- 
lattice structure of the center spins shifts by one site along the leg direction. In order to take 
account of this miss match of the periodicity, we set up three type of the block Hamiltonian. 
Namely, we perform almost the same computation three times and keep the consistency of 
the sublattice structure in the renormalization process by rotating the block Hamiltonians. 

In Fig. 9, we depict the magnetization for A = 1.0. Here, we find a magnetization with a 
plateau at 1/3 of the full moment, which is similar to what we found in the classical model. 
Here we find that the xy-component of the chirality decreases with the field, and disappears 
at the 1/3 plateau. When the field further increases, CSxy appears again. This observation is 
the same as those we found in the section 3. Here we also find the relatively large z-component 
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Fig. 8. Procedures of PWFRG 




(c) (d) 

Fig. 9. The field dependence of (a) the magnetization to, (b) CK^y, (c) CK^, and (d) CSxy as 
functions of the field H for the Heisenberg model A = 1.0. 

of the chirality. This may be due to the fact that we observe the local correlation function 
instead of the order parameter itself. When we reduce A to 0.8 (see in Fig. 10), qualitatively 
similar behavior is observed. It should be noted that, however, the plateau width decreases, 
and the z-component of the chirality increases. For A = 0.7, we find a low field phase with 
z-component of the chirality and it changes discontinuously to a high field phase of the non- 
collinear structure without chirality. (Fig. 11) This is similar to what we found in the classical 
model with a weak XY anisotropy. We also find that the plateau disappears and a jump 
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Fig. 11. The field dependence of (a) the magnetization m, (b) CK^y, (c) CK^-, and (d) CSxy as 
functions of the field H for the Heisenberg model A = 0.7. 



12/15 



J. Phys. Soc. Jpn. 



Full Paper 




(c) (d) 

Fig. 12. The field dependence of (a) the magnetization m, (b) CK^y, (c) CK^, and (d) CS^y as 
functions of the field H for the Heisenberg model A = 0.6. 

appears instead of the plateau. (Fig.ll(a)) The magnitude of the anisotropy at which the 
plateau disappears obtained in this study is close to that obtained by A. Honecher et al?^ {A 
= 0.76 ± 0.03), which is the result of the exact diagonalization study up to 36 sites. As the 
anisotropy becomes stronger, the umbrella structure is stabilized even in higher field region. 
For large anisotropy {A = 0.6), we have no phase transition as the classical case. (Fig. 12) 

Summarizing the dependence of the field-induced phases on the degree of anisotropy, we 
illustrate in Fig. 13 a phase diagram in the {A, H) coordinate. In this figure, we find very 
similar structure to that found in the classical model.^ The boundary between the umbrella 
phase to the Y-shape phase is not determined. For this determination we have to study which 
of Kz or nxy dominates by studying the distance dependence of the correlation functions, 
although so far we obtained only correlation functions between the triangles denoted by crosses 
in Fig. 7. 

6. Summary and Discussion 

From the present observations we found that the ground state of the S = 1/2 XY-like 
Heisenberg antiferromagnetic model depends on the anisotropy sensitively, where the ground 
state configuration is selected by the quantum correction. Then we found that the effect of 
quantum fluctuations in the ground state is very similar to that of thermal fluctuations at 
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Fig. 13. The ground state phase diagram in the {A, H) coordinate. The solid circles represent the 
boundary of the umbrella phase and non-coUinear phase. The triangles represent the upper and 
lower edges of the 1 /3 plateau. 

finite temperatures. In fact, the result for the three-leg ladder shows that the structure of the 
phase diagram (Fig. 13) is very similar to that of the corresponding classical model. In the 
case of two-leg ladder, however, quantum fluctuations cause a special property which has no 
classical counter part. This kind of change of the effect of quantum fluctuations with the width 
of the lattice is also an interesting observation. We note that the difference of the quantum 
and thermal fluctuations is an important future problem. 

Recently, a material with spin-tube consisting of alternate triangles has been found ; it has 
the same lattice as the 3-leg ladder lattice studied in the present paper (Fig.7).^^ We expect 
that the properties found in the present paper, such as the successive phase transitions and 
also their dependence on the spin anisotropy, will be observed in the variety of such materials 
experimentally. 

This work is supported by the Grant-in-Aid from the Ministry of Education, Culture, 
Sports, Science and Technology, and also by NAREGI Nanoscience Project, Ministry of Edu- 
cation, Culture, Sports, Science and Technology, Japan. One of the authors(K.O) is supported 
by Grant for Promotion of Niigata University Research Projects. The simulations have been 
carried out by using the computational facility of the Super Computer Center of Institute for 
Solid State Physics, University of Tokyo. 
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